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The depinning of a contact line is studied as a dynamical 
critical phenomenon by a functional renormalization group 
technique. In D = 2 — e interface dimensions, the rough- 
ness exponent is ^ = e/3 to all orders in perturbation theory. 
Thus, C, = 1/3 for the contact line, equal to the Imry-Ma es- 
timate of Huse for the equilibrium roughness. The dynamical 
exponent is z = 1 — 2e/9 -I- O(e^) < 1, resulting in unusual 
dynamical properties. In particular, a characteristic distor- 
tion length of the contact line depinning from a strong defect 
is predicted to initially increase faster than linearly in time. 
Some experiments are suggested to probe such dynamics. 



Wetting phenomena and contact lines (CLs) appear in 
many manufacturing processes that involve the spread- 
ing of a hquid on a sohd surface R]. Some degree of 
control over the spreading of the liquid and the corre- 
sponding creep of the CL is needed to optimize the de- 
sired characteristics of such processes. In particular, it 
is important to know the effect of surface roughness and 
contaminants at microscopic to mesoscopic scales on CL 
dynamics. Here, we study such effects at length scales 
from 10~^ down to 10~^ centimeters. The upper length 
scale is set by the droplet size or the capillary length (due 
to gravitation) , while the lower length scale is determined 
by the characteristic size of the microscopic defects. 

Surface impurities lead to CL hysteresis^ i.e. a fi- 
nite force is needed to start the fluid spreading. Re- 
cently, the scaling exponents of a driven elastic inter- 
face subject to quenched impurities near a similar depin- 
ning treshold have been calculated through a functional 
renormalization-group (RG) treatment close to four in- 
terface dimensions 0,0]. Here we apply this method to 
calculate various scaling exponents for the slowly advanc- 
ing contact line. The distinction between the two cases 
is that the CL is the termination of the liquid-vapor in- 
terface. We shall assume that the partially wetting fluid 
spreads sufficiently slowly on a heterogenous surface that 
the liquid- vapor interface evolves adiabatically, i.e., the it 
responds to changes in the CL shape instantenously. In 
this case, fluctuations of the CL around its time-averaged 
value reflect the competition between impurities on the 
solid surface and the liquid-vapor surface tension. 

Consider a wetting front on a heterogenous surface in 
the x-y plane with the average orientation of the contact 
line in the x direction as shown in Fig. . In equilibrium 
on a homogenous interface, the macroscopic contact angle 



is determined by the Young condition, 
isv - ISL - 7 COS 9 = 0. 



(1) 



In the above, 7sy, Jsl, and 7 are the interfacial ten- 
sions for the solid-vapor, solid-liquid, and liquid-vapor 
interfaces, respectively. The heterogeneities (aka defects) 
on the surface can be modeled as fluctuations in the dif- 
ference of local interfacial energy densities, 

f{x, y) = -isv{x, y) - jsl{x, y) - {isv - Isl)- (2) 

The average of / is zero, while its correlations satisfy. 



(/(x,2/)/(x',2;')) = A(r/a), 



(3) 



where r — {x — x') + (jj — y')\ a is the characteristic 
size of defects, and A is a function that decays rapidly 
for large values of its argument. (For a self-afRnely rough 
surface, A may have a slow algebraic decay [Q. Such 
situations will not be explored here.) Fluctuations of the 
CL from its average position, hcL{t) = vt, are denoted 
by h{x,t), and thus satisfy 



{h{x,t)) =0. 



(4) 
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FIG. 1. Geometry of the system. 



(The overline denotes a time average.) The capillary 
energy associated with small deformations of a CL was 
calculated by Joanny and de Gennes |j] (in the 8-^0 
limit) as 
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where h{q) is the fourier transform of the contact line 
profile at a given time, and a and L are the lower and 
upper cutoff length scales mentioned earlier. The unusual 
q-dependence of the energy functional, and the resulting 
nonlocal dynamics, reflects the fact that perturbations 
of wavevector q on the CL induce deformations into the 
liquid-vapor interface up to a distance |<z|~^. 

A random contribution to the liquid-surface energy 
comes from the defects fllf , 



Ur, 



dx 



hcL+h{x) 



dy f{x,y). 



(6) 



If the heterogeneity is strong enough, there is contact 
angle hysteresis ||l|,Q . This arises from many metastable 
configurations of the CL profile, which are given by local 
minima of the free energy. The surface tensions 7sl , Isv 
in Eq.(yj) must then be interpreted as spatial averages 
over the CL position. Thus, the macroscopic contact an- 
gle will depend on the particular CL profile. If a macro- 
scopic force of F per unit length is exerted on the CL, 
it will move with a finite velocity v only for F > Fc- 
This occurs when the metastable state with the largest 
8, usually referred to as the advancing angle 8a [|l],||, 
becomes unstable. For small capillary numbers, 8 ~ 8a. 
As the CL slowly advances, energy is dissipated 
through various mechanisms ||l| . If the dissipation in the 
vicinity of the CL dominates for small v, the drag force at 
a point x on the CL is simply related to the local velocity 
V + dth(x,t) through a microscopic mobility fj,. In this 
case, the equation of motion for the CL is obtained by 
equating the drag force to the applied force, —6U/6h{x), 
as 



dh{x,t) 
dt 
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f[x,vt + h{x,t)] +F. 
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In the above formula, v is determined self-consistenly by 
enforcing Eq. (0). Note that Eq.(R) is only valid for 
w > 0, and F ^ F^ from above when v -^ 0. The 
dynamical contact angle 8 is given by the force balance 
equation, 



F ~ Fc — 7(cos8a — cos 8). 



(8) 



A recent experiment on depinning from a single defect 
pI is in agreement with the deterministic (/ = 0) form 



of this dynamics: Upon depinning from the defect, the 
time evolution of the CL profile is given by 



h{x, t) '^ In 



(ct? 



L2 



(9) 



where the width of the logarithmic profile increases lin- 
early with a characteristic velocity c = jQ'^fj, after depin- 
ning. 

By analogy to similar previously studied systems 
|g|,||J3] , for F slightly above treshold, we expect the aver- 
age velocity of the contact line to scale as 



AiF - Fc) 



13 



(10) 



where (3 is the velocity exponent, and A is a nonuni- 
versal constant. Superposed on the steady advance of 
the CL are rapid "jumps" as portions of the line depin 
from strong pinning centers. Such jumps are similar to 
avalanches in other slowly forced systems, and have a 
power-law distribution in size, cut-off at a characteristic 
correlation length ^. On approaching the threshold, ^ 
diverges as 



e ^ (F - F, 



(11) 



defining a correlation length exponent v. At length scales 
smaller than ^, the interface is self-affine, i.e. with corre- 
lations satisfying the dynamic scaling form 

((Mx, t) - h{x\ t')f) ^{x- xTg (^ ^^Zty ) ' (12) 

where C and z are the roughness and dynamic exponents, 
respectively. The scaling function g goes to a constant as 
its argument approaches 0, C is the wandering exponent 
of an instantenous CL profile, and z relates the average 
lifetime of an avalanche to its size by r(^) ~ ^^. 

In many aspects, Eq.(M) is similar to the model ana- 
lyzed by Narayan and Fisher using the formalism of Mar- 
tin, Siggia, and Rose w] (MSR), through an expansion 
around a mean field solution 0]. To use this method, 
it is necessary to generalize to a D-dimensional interface 
for a systematic computation of the exponents through 
an e expansion. Introducing an auxiliary field h{x, t), the 
MSR generating functional is, 

Z = [dh] [dh] expli / d^x dt h{x, t)T[h, /] I , (13) 
T [h, f] = A^-i (at/i(x, t)+v)- /[x, vt + /i(x, t)] - F 
+78^ f d^x'KDix-yi')h{x',t). 

In the above, Koiy) (^ |x|"(^+i) for a < |x| < L) is the 
generalized interaction kernel in D dimensions. Mean 
field (MF) theory is obtained by replacing the capillary 
forces on a portion of the CL with Hookean springs. This 
gives 



M 



^(dhM^^\^_ 



V dt 



-khMF{t) + f[vt + hMF{t)] + F, 

(14) 



independently for each x, and identical to the mean field 
model analyzed by NF [||j^. Here again, v{F) is deter- 
mined by the self-consistency condition {Hmf) — 0. 

The MF solution depends on the type of irregularity: 
For smoothly varying f{x,y), j3mf = 3/2, whereas for 
"mesa" defects, i.e. cusped f{x,y), Pmf = 1 Q^^- Some 
recent experiments on reasonably clean surfaces by Strom 
et al. [g| are in agreement with Eq.([lO|) for /3 = 3/2 in 
the low capillary number regime. Even though it is gen- 
erally believed that the scaling relation in Eq.(]lO|) holds 
|10| , it is clear that the prefactor is nonuniversal, and 
various experiments to date have obtained widely vary- 
ing v-F behaviors |11| . To our knowledge, there are no 
experiments so far that have systematically investigated 
strongly heterogenous surfaces with reasonable accuracy. 

Following the treatment of NF am , we use the mean 
field solution for cusped potentials, anticipating jumps 
with velocity of 0(1), in which case (3m f = 1- Af- 
ter rescaling and averaging over impurity configurations, 
we arrive at a generating functional whose low-frequency 
form is 



Z= l[dH][dH]exp{-~S), 



S = / d^xdt [F - Fmf{v)] i?(x, t) 



duj d^q 
2^ {2n)D 



H{-q,-uj){-iuj+\q\)H{(i,uj) 



— / d°y.dtdt' H{K,t)H{-K,t') 

xC[vt~vt' + H{x,t)-H{yi,t')]. (15) 

In the above expressions, H and H are coarse-grained 
forms of h and h, in the sense that response and corre- 
lation functions of h, h are the same as those of H, H 
at low frequencies. F is adjusted to satisfy the condition 
{H) = 0. The function C{vt) is initially the mean-field 
correlation function {hMF{t)hMF{t + r)). 

All non-gaussian terms in the generating functional Z, 
including ones that are not shown in Eq. (|l^) , decay away 
at large length and time scales for D > D^ = 2. Thus, 
for D > Dc, exponents can be easily determined from 
the gaussian theory as zq = ^o = /3o = 1- The interface 
is flat at long length scales {C,o = (2 — D)/2). 

At D — Dc, there are an infinite number of marginal 
operators that can be cast into a single marginal function 
for 11 — > 0. The RG is carried out by integrating over a 
momentum shell and all frequencies, followed by a scale 
transformation x — > bx, t — > b^t, H -^ b'^H, and H — > 
b^^^H, where b = eK For D < Dc, there are a number 
of exact exponent relations M : 



iy = i/{z + e), 
^ = 1/(1-0- 



(17) 
(18) 



The first result comes from the scaling of nonlinear re- 
sponse functions, requiring H to scale as vt. The second 
relation follows from the renormalization of i^ — F]\ip{v) 
through the first term in Eq.(15). Finally, the last re- 



lation results from an exact statistical symmetry of the 
system that fixes the form of the static (w = 0) response 
function {dh{x')/df{x)). These relations determine all 
exponents in terms of C, and z. The diagrams appearing 
in the renormalization of C{u) are exactly the same as in 
Ref . H] , with the exception of the form of the bare prop- 
agator, which is (— iw -I- |g|)~"^ instead of {—iu + q^)"^ . 
Thus, the renormalization of C in D = 2 — e interface 
dimensions, computed to one-loop order, gives the recur- 
sion relation, 



dC{u) 
dl 



[e + 2e + 2{z-l)]C{u) + (:u 
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du 
dCju) 
du 



(19) 



NF showed that all higher order diagrams contribute to 
the renormalization of C as total derivatives with respect 
to u, thus, integrating Eq.(|l9|) at the fixed-point solution 



dC*/dl = 0, together with Eqs.(0) and i 
e/3 to all orders in e, provided that / C* 
case of a CL (D = 1), we thus obtain 



C = 1/3, and 



3/2. 



y), gives C = 
^ 0. For the 



(20) 



This value of the roughness exponent coincides with the 
Imry-Ma estimate of Huse given in Ref. 0| for the equi- 
librium roughness. This is a consequence of the fact that 
C{u) remains short-ranged upon renormalization, imply- 
ing the absence of anomalous contributions to C- 

To calculate the dynamical exponent z, we need the 
fixed-point function C* , which is obtained only to 0(e) 
from the above analysis. Furthermore, the calculation 
requires some knowledge of the high-frequency form of 
response and correlation functions, i.e., a low- frequency 
analysis of Eq. ( p^ is insufficient to describe scaling prop- 
erties of the system. Nevertheless, NF calculate the ex- 
ponent 9 to 0(e) as 9 + C = {e- C)/3 + O^e^). The 
computation for the CL gives the same result, yielding 



z = l-2e/9 + 0(e2), 
13= l-2e/9 + 0(e2). 



(21) 
(22) 



P={z- Ov, 



(16) 



Note that (3 < z even though the difference is O(e^). Nat- 
termann et al |^ obtain the same results to 0(e) by di- 
rectly averaging the MSR generating function in Eq. (Oh . 
The treatment of NF has the advantage of expanding 
around a solution with Fc ^ 0, resulting in a better be- 
haved theory. In particular, systems that are described 
by vector fields instead of scalar fields, like the treshold 
critical dynamics of flux lines in three dimensions, require 



the use of their approach, since the expansion has to be 
made around a solution with Fc ^ Q . 

The above scahng exponents describe a CL advancing 
at low capillary number on a dirty surface, and also apply 
to surfaces with microscopic roughness, i.e. with short 
range slope-slope correlations. The roughness exponent 
of C = 1/3, equal to its equilibrium value, can be directly 
measured by microscopic imaging of a slowly advancing 
CL. A number of macroscopic experimental tests are also 
possible: the velocity of the CL is given by by u ~ (F — 
Fc)^ with /3 < 1 (See Fig. ). The root mean square width 
of the CL profile should increase like W ^ ^'' ^ y^C/i^P 
as the threshold is approached. These relations can be 
checked through tensiometric measurements ||ll[, where 
the capillary force on a solid immersed into a liquid is 
measured directly. 



XT = 



d{Jh) 




FIG. 2. Predicted velocity-force behavior for a 
depinning contact line; (3 =7/9 to first order in e 



=2-D. 



A dynamical exponent of z < 1 suggests that the relax- 
ation of the CL is very different on a dirty surface. Upon 
depinning from a strong defect, the width of the loga- 
rithmic CL profile initially grows faster than linearly in 
time, in other words, with a characteristic velocity that 
increases with time: c{t) ^ t^^^^^l^ . This is, of course, 
not physical at arbitrarily large length and time scales. 
We have assumed that the liquid-vapor interface retains 
its equilibrium shape determined by the CL profile at 
all times. This picture will break down as characteristic 
velocities become comparable to capillary wave velocity 
Ccap = \Jl I P of the liquid-vapor interface. The scaling 
regime in between should still be accessible to experi- 
ments in which the depinning from a single strong defect 
is observed on a dirty or rough surface. 

Finally, NF carried their analysis further to speculate 
exponents helow the treshold when the driving force is 
increased monotonically towards i^c 0. In particular, 
they postulated an avalanche distribution 



dF^ 



(Fc-F)- 



(24) 



If i/_ = v, and the scaling hypothesis holds for the CL, it 
then follows that k — D — 1/v^, and 7=1-1- C,v-. These 
exponents are then given exactly by 7 = i^_ = 3/2, and 
K = 1/3 in D = 1. 
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Prob(width of avalanche > £) « —p{£/^-), (23) 

£k 



with ^_ ~ {Fc — Fy- , and the mean polarizability 



